Coherence Properties of Guided- Atom Interferometers 



m 
o 

O 
(N 

Oh' 
D 



O 



I 

O 

o 



(N 
(N 
0^ 
O 
m 
o 



o 



H. Kreutzmann,^'Q U. V. Poulsen/ M. Lewenstein,^ R. Dumke,^ W. Ertmer,^ G. Birkl,^ and A. Sanpera^ 

''institut fiir Theoretische Physik, Universitat Hannover, 30167 Hannover, Germany 
^Institut fiir Quantenoptik, Universitat Hannover, 30167 Hannover, Germany 
(Dated: February 2, 2008) 

We present a detailed investigation of the coherence properties of beam sphtters and Mach- 
Zehnder interferometers for guided atoms. It is demonstrated that such a setup permits coherent 
wave packet splitting and leads to the appearance of interference fringes. We study single-mode and 
thermal input states and show that even for thermal input states interference fringes can be clearly 
observed, thus demonstrating the multimode operation and the robustness of the interferometer. 
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The investigation and exploitation of the wave proper- 
ties of atomic matter is of great interest for fundamental 
as well as applied research and constitutes, therefore, one 
of the most active research areas in atomic physics and 
quantum optics. Of special interest is the field of atom in- 
terferometry 1] . For an interferometer, it is crucial that 
the beam splitters and mirrors are coherent, i.e., they 
must not disturb the phase of the matter wave in an un- 
controllable fashion. Then a phase shift in one of the 
paths results in a change of the output signal, and any 
external influence inducing a phase shift is, in principle, 
accessible to measurement. Compared to light interfer- 
ometry, matter wave interferometry with cold atoms of- 
fers a much higher sensitivity for certain applications . 
Furthermore, atoms couple efficiently to a wider variety 
of external interactions, thus extending the applicability 
of interferometric measurements 0. 

A novel approach arises from the use of guided 
atoms 0) S H- Miniaturized setups for matter wave 
interferometry with increased stability, large beam se 
aration, and large enclosed areas become possible [" 
0, . These features are specifically appealing to atom- 
interferometrical measurement of inertial forces Q and 
to the investigation of Bose-Einstein condensates in mi- 
crostructures |2| . Due to the physics involved in the guid- 
ing and beam splitting processes, the construction and 
analysis of guided-atom interferometers become challeng- 
ing tasks. In order to assess the level of performance that 
can be reached with realistic setups, coherence and inter- 
ference, also for mixed input states as well as non-perfect 
beam splitters, have to be investigated in detail. 

In this Letter we study guided-atom interferometers 
for neutral atoms by solving numerically the time depen- 
dent Schrodinger equation for realistic, experimentally 
accessible configurations. Our study addresses the main 
issues of atom waveguide propagation, as well as coher- 
ence and interference using X-shaped guided-atom beam 
splitters 0, 0,0, 113 • particular, our calculations a pply 
to the interferometer scheme originally proposed in [llj 
in which neutral atoms are guided in dipole potentials 
created by micro-fabricated optical systems. Many of 
the experimental prerequisites for that proposal such as 



atom guides, beam splitters and even geometries compos- 
ing a complete Mach-Zehnder interferometer for atoms 
have been realized already 

In the scheme considered here, the beam splitters are 
achieved by crossing two optical waveguides at some an- 
gle 0, H • Such a beam splitter splits the atomic wave 
packet in coordinate space without affecting the internal 
state. Nevertheless, internal state selective interferome- 
try is also possible in this scheme |3| . During the splitting 
process the system might exhibit quantum reflections and 
tunneling between adjacent guides and, therefore, the 
dynamics is in general complicated. We will, however, 
demonstrate that such a beam splitter is coherent even 
for a thermal distribution of atoms with an average en- 
ergy far exceeding the level spacing of the transverse con- 
finement. 

To study the properties of the beam splitter and the 
interferometer we use the Split- Operator method to solve 
the time dependent Schrodinger equation 
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where the potential V{x,y) includes the waveguides po- 
tentials as well as any other relevant potential in the 
problem. For simplicity we shall assume throughout this 
paper that the atomic wave packet is tightly confined 
in the third dimension so that the dynamics is well de- 
scribed within a two-dimensional treatment. 

The central element of any interferometer is the beam 
splitter. Consequently, we start our discussion with a de- 
tailed analysis of an X-shaped guided-atom beam splitter 
(BS) created by crossing two identical waveguides and 
hj at an angle 7: 

UBs{x,y) = C/i(a;, j/)-|-t/,(a;cos7-?/sin7,?/cos7-|-a;sin7) 

Each waveguide consists of a Gaussian potential of depth 
Uq and width a; thus a waveguide along the y-direction 
centered at a;o is represented by: 
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Coupling from one guide to the other occurs for any an- 
gle 7 ^ 90°. We observe, however, that the doubling of 
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the potential depth at the crossing of the two waveguides 
induces quantum reflections and a highly non-adiabatic 
dynamics for typical initial momenta of the atomic wave 
packet. These effects are clearly undesirable for interfer- 
ometry. For a micro-optical realization of the waveguides, 
they can be easily avoided by reducing the light intensity 
at the crossing through overlay of an absorptive mask or 
by adding a compensating extra potential (e.g., a blue 
detuned laser field) so that the depth of the total poten- 
tial at the crossing is equal to the one of each waveguide 
alone. In our simulations we have taken into account that 
this compensation might not be perfect. 

Our choice of parameters closely matches the rele- 
vant experimental parameters of refs. H 113 for 85 Rb 
atoms guided in dipole potentials, far detuned below the 
5Si/2 — > 5P3/2 transition at 780 nm. A typical configu- 
ration consists of waveguides of width a — 0.54 /im (cor- 
responding to a Gaussian beam with waist wq = 1.1 
/im) at a laser wavelength of 830 nm and an intensity of 
/ = 1.1 X 10^ W/cm^ (less than 1 W of required laser 
power). The depth of potential is Uq = 75 /iK, the ground 
state vibrational frequency oj — 160 x 10'^ s~^, and the 
rate of spontaneous scattering Fs = 2.6 (thus it can 
be neglected in our discussion). In our simulations, as 
initial state we consider an atomic wave packet located 
in one of the waveguides (Li in Fig. ^ at a typical dis- 
tance of 2.5 /im from BSi. We perform simulations for 
both single and multimode transverse occupation of the 
waveguide. In general, the initial transversal state can 
be described by a thermal mixture 

where | n) denotes the n**^ eigenstate of the waveguide po- 
tential, En denotes its energy and Z ensures proper nor- 
malization. In the experiment Q an ultra-cold atomic 
sample with temperature of T = 20 /iK is loaded into 
the waveguide from a single dipole trap. Our set of pa- 
rameters results in a mean transverse occupation num- 
ber of ( n ) ~ 16 in the waveguide. For both single-mode 
(p = I I) and multimode cases, we assume that the 
initial atomic wave packet has a Gaussian profile along 
the longitudinal direction of the waveguide with a mean 
momentum py and a spread of Apy . Values of the mean 
momentum are in the range of 5 — 10 recoil momenta Pr 
(pr = \J1hmLij-[^ Wr — 2A X 10'^ s~^), with a momentum 
spread Apy = 2pr. 

The splitting of the initial wave packet depends on its 
initial transversal states, its initial longitudinal momen- 
tum Py, and on the angle 7 between the guides. In order 
to achieve efficient deflection in the beam splitter, the 
atomic wave packet has to spend enough time at the in- 
tersection. Defining the crossing time as tc = cr/vy with 
Vy ~ py/m, efficient splitting requires tc ^ Ti/En- For a 
fixed angle 7 and a fixed initial momentum pj,, we observe 




FIG. 1: Schematic view of a Mach-Zehnder interferometer 
for guided atoms using 4 identical waveguides crossing at an 
angle 7. The initial atomic wave packet is represented by a 
dot in waveguide Li below the first beam splitter (BSi). The 
location of the phase shift potential U{x,y) is also depicted. 

that the dependence of the splitting ratio on the initial 
transverse state | n) is very strong. This is displayed in 
Fig. 13 where we plot the fraction of atoms transmitted 
towards BS2 (T) and the fraction deflected towards BS3 
(D) as a function of the initial transverse state for 7 = 45° 
and two different initial longitudinal momenta, py = lOpr 
andpj, = 5pr. To "count" the number of deflected (trans- 
mitted) atoms we use an absorbing box at L3 (Li) after 
the beam splitter and integrate the loss of norm in each 
box with time. The transmitted and deflected fraction 
do not always add up to unity. The missing fraction con- 
sists of atoms backscattered into the downward sections 
of waveguides Li and L3. For py — lOpr, an approxi- 
mately 50/50 splitting ratio occurs for transverse initial 
states with quantum numbers n ~ 8 — 11. Losses due to 
backscattering are very small in this case. For py = 5pi. 
the optimal splitting ratio occurs for n ~ 2 — 3 although 
losses are now significantly higher. In both cases, the de- 
flected fraction is narrowly peaked around its maximum 
evidencing that the beam splitter acts as filter for trans- 
verse states. This implies that even for thermal input 
states, efficient splitting occurs only for a narrow group 
of states around the optimal one. In both cases the frac- 
tion of deflected atoms is very small for the ground state 
(p= |0)(0|) of the potential. 

The following simple picture helps to understand the 
selection of the optimal state | n) for a given initial mo- 
mentum Py and angle 7: At the intersection, the wave 
packet is no longer conflned transversally and, therefore, 
expands according to its initial transverse momentum 
distribution. Optimal splitting occurs for typical trans- 
verse momenta px = ^J2mEn fulfilling px — Pyt&Ti^. 
Approximating the center of the waveguide by a har- 
monic potential of frequency uj one obtains the estimate 
Tiopt — Py tan^(7)/(2m7iw). Lower transverse states can 
be made optimal by lowering py and/or lowering 7. A 
lower rtopt also implies fewer states within the peak of 
non-negligible splitting ratios thus enhancing the filter- 
ing effect of the beam splitter. Despite the simplicity 
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FIG. 2: Splitting efficiency of guided-atom beam splitter: 
transmitted (T) (circles) and deflected (D) (triangles) frac- 
tion of the incoming wave packet as function of the initial 
transverse mode of the waveguide for 7 = 45°. Full symbols 
correspond to py = lOpr and open symbols to = 5pr- 
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FIG. 3: Interference fringes for a single-mode (p — |0){0|) 
initial state versus the strength d of U{x,y). (a) Fraction of 
atoms (in % of the input atom number) in the top output, 
(b) Fraction of atoms (in %) in the top output for the ini- 
tial interferometer size (circles) and the optimized size (L2 
displaced in x-direction) (stars). 



of the argument, it agrees well with our results depicted 
in Fig. 121 and with our simulations for different initial 
longitudinal momenta Py and angles 7. Only for angles 
7 < 25° we find strong deviations from our estimate since 
there backscattering and tunneling start to play an im- 
portant role in the splitting dynamics and the simplified 
argument fails to reproduce the numerical results. 

The coupling between transverse and longitudinal de- 
grees of freedom in the splitting process results in a com- 
plex dynamics after the beam splitter. We calculate the 
longitudinal phase profile of each wave packet before it 
reaches the next beam splitter. The transmitted wave 
packet propagates along the waveguide Li almost with- 
out distortion and its phase profile agrees well with the 
phase profile of a wave packet propagating freely in a 
transverse harmonic potential 12j. The motion along L3 
is more involved. Although the overall behavior of the 
phase resembles the one of the wave packet along Li, 
even for a single mode initial state, the wave packet ex- 
hibits transverse oscillatory motion of frequency lu. This 
is caused by the excitation of a coherent superposition of 
transverse modes and clearly demonstrates that the split- 
ting process is non-adiabatic. This oscillatory motion has 
to be taken into account for optimizing the efficiency of 
the full Mach-Zehnder interferometer. 

The full interferometer is realized by appropriately 
crossing four identical waveguides. Due to numerical lim- 
itations, we constrain the distance between Li and L2 to 
about 7 ^m. |0| Atoms loaded in the waveguide Li 
with initial momentum py propagate to the first beam 
splitter BSi and are then guided through Li and L3 to- 
wards BS2 and BS3 respectively. Here BS2 and BS3 act 
as mirrors, but due to the additional output ports also es- 
tablish new loss channels. Note that this Mach-Zehnder 
configuration is not symmetric even for 50/50 beam split- 
ters, since, as we have already discussed, the dynamics in 
each arm of the interferometer becomes quite distinct af- 
ter the first beam splitter. Following splitting processes 



in the two arms are no longer equivalent. This lack of 
symmetry excludes the possibility of achieving 100% vis- 
ibility (defined as the amplitude [i.e., 1/2 (max-min)] of 
the modulation in the top output divided by the average 
in the top output) and demands optimization strategies 
to improve the visibility. 

To demonstrate that coherence is preserved during the 
propagation through the full interferometer, we calcu- 
late the output signals after the final beam splitter BS4 
as a function of the depth d of an additional potential 
U{x,y) = _d[/oe-(^-*«)'/2'^'e-(''-*'»)'/2<T2 which is in- 
serted in one of the arms of the interferometer (between 
BS3 and BS4) to induce a phase shift between both arms. 
This extra potential, which can be either attractive or 
repulsive depending on the sign of d, has a Gaussian pro- 
file along the waveguide and smoothly lowers {d > 0) 
or increases {d < 0) the potential depth of that part of 
the waveguide. Our results are summarized in Fig. |21 
and Fig. ^ Clear periodic modulations in the number 
of atoms exiting each output port (labeled top and side 
outputs for clarity) appear as a function of d for both 
single-mode (Fig.OJ and thermal initial state (Fig. 0)). 

For the single mode case, we choose, for simplicity, 
the ground state with longitudinal momentum py = lOp^ 
as initial state. In Fig. Ol^a), we display the fraction of 
atoms at the top output versus d. The expected modu- 
lation of the atom number as a function of d is clearly 
visible. We compare the period of the oscillation with 
the one obtained using a simplified model of the phase 
shift (j) introduced by the additional potential U using the 
classical action S — JdtC, where C is the Lagrangian. 
Despite the oversimplification of this model, the calcu- 
lated period is only a few % too small for low values of d. 
For larger d the increase in the oscillation period is more 
accurately described using a WKB approximation. 

The combined atom number of both outputs of the 
interferometer is, for this case, small (sa 1% of the ini- 
tial atom number) since, as shown in Fig. |21 the ground 
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FIG. 4: Interference for single-mode and thermal input states: 
fraction of atoms (in %) in the side (triangles) and top output 
(circles) versus the strength d of U{x,y). Filled symbols cor- 
respond to the transverse state with optimal splitting n = 2, 
open symbols to an initial thermal state at T=20 /iK. 



state splits very inefficiently and most of the atoms leave 
the interferometer at BS2 and BS3. The visibility of 
the fringes is also low 2%). Better visibility can be 
achieved by choosing as the initial state the one with the 
optimal splitting ratio (c.f. Fig. 0)) and using an opti- 
mized geometry. An optimized geometry aims at equal- 
izing as much as possible the contributions from both 
arms of the interferometer to the final signal. This can 
be accomplished by choosing the distance between the 
waveguides Li and L2 such that the wave packet that 
propagates in L3 with transverse oscillations enters the 
beam splitter BS3 with a mean transverse momentum 
appropriate to maximize the fraction of atoms deflected 
towards BS4. This optimization immediately leads to a 
visibility of ~ 7% as shown in Figl^Ib). 

Experimentally, a single transverse mode wave packet 
is difficult to achieve. A more realistic scenario corre- 
sponds to an initial thermal occupation of the transverse 
modes (Eq.O. For this case, we calculate the final output 
signal as a classical (Boltzmann) weighted average of the 
signal obtained for each transverse mode | n) separately. 
In our calculations, we now use an optimized geometry 
but we consider, nevertheless, that the compensation of 
the double potential in the beam splitters is not perfect 
allowing for a 5% mismatch in potential depth at the 
beam splitter. In Fig. 0] we display the output signal for 
an initial state corresponding to all atoms in the optimal 
spliting state p = | 2) (2 | (c.f. Fig. and for a thermal 
(multimode) state corresponding to T=20 fiK . In both 
cases py = 5pr and 7 = 45°. Since the beam splitter 
is state selective, interference fringes are mostly due to 
states with a good splitting ratio as evidenced by the fact 
that maxima and minima appear approximately at the 
same positions in both cases. The combined atom num- 
ber for the optimal state is about 24% of the initial atom 
number, with visibilities up to 15%. (Higher visibilities 



of « 23% can be obtained for p = | 1) (1 1 but with a much 
lower signal). For the thermal initial state the different 
splitting ratios corresponding to the different transverse 
states lower the total output signal to « 17% but with 
visibilities up to 10%. This clearly demonstrates the per- 
sistence of coherence and interference even for thermal 
input states and allows for the operation of the inter- 
ferometer as a multimode device. In fact, compared to 
the zero temperature case (ground state), a thermal state 
which inherently contains the optimal state will dramat- 
ically improve the performance of the interferometer for 
large 7. 

In summary, we have shown that guided-atom con- 
figurations allow for significant extensions of matter wave 
interferometry. An X-shaped beam splitter created by 
crossing two identical waveguides (with the doubling of 
the potential properly compensated) preserves coherence 
and acts analogously to a "color filter" , selecting only a 
few optimal radial states which contribute to the final 
signal. In this way coherence is preserved throughout 
the interferometer. As a consequence, temperature is not 
necessarily a limitation for the observation of interference 
fringes but can rather be a requisite to ensure a good 
performance. Our simulations show that an additional 
potential of variable strength in one of the arms of the 
interferometer can produce a straightforward proof for 
coherence and interference. Finally, we have also demon- 
strated improvements by applying optimization strate- 
gies. 
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